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ABSTRACT 

We present a structural analysis of NGC 891, an edge-on galaxy that has long been 
considered to be an analogue of the Milky Way. Using starcounts derived from deep 
HST/ACS images, we detect the presence of a thick disk component in this galaxy with 
vertical scale height hz = 1-44 ± 0.03 kpc and radial scale length Iir — 4.8 ± 0.1 kpc, 
only slightly longer than that of the thin disk. A stellar spheroid with a de Vaucouleurs- 
like profile is detected from a radial distance of r ~ 0.5 kpc to the edge of the survey 
at r ~ 25 kpc; the structure appears to become more flattened with distance, reaching 
q = 0.50 in the outermost halo region probed. The halo inside of r ~ 15 kpc is 
moderately metal-rich (median [Fe/H] ~ —1.1) and approximately uniform in median 
metallicity. Beyond that distance a modest chemical gradient is detected, with the 
median reaching [Fe/H] ~ — 1.3 at r ~ 20 kpc. We find evidence for subtle, but very 
significant, small-scale variations in the median colour and density over the halo survey 
area. We argue that the colour variations are unlikely to be due to internal extinction 
or foreground extinction, and reflect instead variations in the stellar metallicity. Their 
presence suggests a startling conclusion: that the halo of this galaxy is composed of 
a large number of incompletely-mixed sub-populations, testifying to its origin in a 
deluge of small accretions. 



1 INTRODUCTION 

Remarkable progress has been made in recent years in un- 
derstanding galaxy formation and evolution. High redshift 
observations have revealed the star formation history of the 
universe and the evolution of gala xy morphology (see, e.g., 
iBell et alj|2005l : iRvan et alJ lioOSV However, the nature of 
look-back observations does not allow one to study the evo- 
lution of individual galaxies, and low-mass or small-scale 
structures remain out of reach in all but the nearest galax- 
ies. Hence observations in nearby galaxies are required to 
complement the samples at cosmological distances to answer 
such fundamental questions as: which building blocks of high 
redshift galaxies end up in what type of local galaxy?; how 
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much of the stellar content of different galactic morphologi- 
cal components is formed in situ and how much is accreted? 

In recent years, it has been increasingly recognised that 
many of the clues to the fundamental problem of galaxy for- 
mation are preserved in stars in the outskirts of galaxies, 
offering the unique opportunity to address these questions 
l|Freeman fc Bland-Hawthorni 2002). Hierarchical formation 
models suggest that galaxy outskirts form by accretion of 
minor satellites, predominantly at early epochs when large 
disk galaxies were assembling for the first time. The size, 
metallicity, and amount of sub-structure in the faint out- 
skirts of present-day galaxies are therefore directly related 
to issues such as the small scale properties of the primordial 
power spectrum of density fluct uations and the supp ression 
of star formation in small halos (|Springel et al.ll2005t ). 

The Milky Way and the Andromeda galaxy (M31) have 
played a pivotal role in studi es of the faint outer regions 
of ga laxies over many decades ijFreeman fc Bland-Hawthornl 
120021 ). Indeed, almost all we know about the properties of 
the faint outskirts of galaxies is based on observations of 
these two galaxies. In recent years, we have seen dramatic 
progress in large scale mapping of the Milky Way, with 
large surveys that are measuring (and will measure) pho- 
tometry, kinematics, and chemical abundances of a large 
number of individual stars over a wide volume of the 
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G alaxy (SPSS -iGunn et al.ll2006l ; llvezic et al.ll2008l; RAVE 



■ ISteinmetz et al.l l200d : SDSS-II /SEGUE - iRockosil 120051 ; 
GAIA; Pan-STARRS). It is no exaggeration to say that the 
study of these data sets, to unravel the processes by which 
the Galaxy came into being, will be among the chief astro- 
physical goals for years to come. 

Although these Galactic surveys will significantly ad- 
vance our understanding of the assembly process, it is un- 
likely that a sample of one galaxy will provide the defini- 
tive solution to the nature and origin of stellar popula- 
tions in galaxies in general. A necessary step to fully ex- 
ploit Galactic surveys and thereby establish a comprehen- 
sive picture for spiral galaxy formation is to determine 
whether the Galaxy is suitably typical by studying other 
giant spirals. It is precisely for this reason that the An- 
dromeda galaxy has been surveyed in gr eat detail in recent 
years both with wide-field photometry l|lbata et alj|2001al ; 
Fcrauson ct al ] |2002l;llbata et al 1120071) and massive multi- 



object spectroscopy dlbata et alj2005tlChapman et "al]|2006l ; 
iKalirai et all 12009 : iGuhathakurta et al.ll200rj ). These anal- 
yses show that the Milky Way and M31 are very different, 
with M31 displaying an abun dance of signs of both ancient 
and recent merger activity (|Block et alj 120061 ; llbata et al.l 
l2007h . while the Milky Way appears to have led an unusually 
quiet existence l|Hammer et al.1 120071 ) . Significant effort is 
currently being devoted to interpret the cosmological impli- 
cations of the substructures observed in the halos of nearby 
galaxies (e.g. jFont et al.l [20081 ; Ijohnston et alll2008l ). 

In the currently favored galaxy formation models, ha- 
los are populated by stars which are tidally-s tripped from 
satell ites as they fall into the host potential l|Abadi et al.l 
2006). The structure, chemistry, and ubiquity of stellar ha- 
los are thus exp ected to reflect the details of the galaxy 
assem bly process l|Bullock fc Johns ton 2005; Johnston et al.l 
2008). The realization that hierarchical models of galaxy for- 
mation almost inevitably predict one to two orders of magni- 
tude more d warf galaxies than are actually seen ar ound the 
Milky Way (|Moore et al.lll999l ; iKrypin et all 119991 ) has led 
to suggestions that most of the dwar fs formed in the ear ly 
universe have dissolved into the halo l|Bullock et al.l l2000), 

Observationally however, the nature and the origin of 
galactic stellar halos remain elusive. Excitement about the 
formation of stellar halos has been sparked in recent years 
after the discoveries of disrupted dwarf galaxies and sub- 
structures in star samples in the Galaxy Jl bata ct al., 1 19941 ; 
lYannv et al.ll2000l; llbata et al.ll2003l; iMartin et al.ll2004r) and 
in M 31 |lbata et al.ll2001al ; IZucker et al.l 12004 llbata et al.1 
120071 ). reinforcing the concept that satellite accretion could 
be an important mechanism for building up galaxies. Stud- 
ies of the inner parts of the M31 "halo" initially revealed 
dramatic differences between the halo populations o f that 
galaxy and those of the Galaxy l|Durrell et al.l 12004) . The 
Milky Way halo seems to be populated by old, metal-poor 
stars, while a few fields in the halo of M31 show a large 
population of inter mediate age stars with a much higher 
overall metallicity l|Brown et al.1 12003. 2007). Red horizon- 
tal branch stars are also found as distant as 40 kpc from the 
nucleus of M31 (|Rich et al.ll2004l) suggesting that the stellar 
halo of M31 underwent considerably more enrichment than 
did the Galaxy halo. However, the conclusion that the M31 
halo is globally younger and more metal-rich than that of 
the Milky Way cannot be drawn from extant data, since 




Figure 1. The locations of the three ACS fields HI, H2 and 
H3, superimposed on photographic data from the Digitized Sky 
Survey. 



the fields in M31 in which the above results were obtained 
have been found to be significantly co ntaminated by vari- 
ous accretion events l|lbata et al.ll2007h . Furthermore, spec- 
troscopy of stars either selected by their kinematics to be 
halo members or sufficiently in the outskirts of the galaxy re- 
veals the presence of a pressure-supported metal-poor popu- 
lation _sjmikMl»tJiepopji^ the Galactic 
halo (|Chapman et al'lbOOrMKalirai et al.ll2006l ). The current 
observational evidence therefore demonstrates that halos are 
complex structures; and larger samples will be required to 
fully understand how they form and evolve. 

To answer the fundamental question of how typical 
are the Local Group massive galaxies, and to constrain 
the assembly history of spirals, we need to probe beyond 
the Local Group, and we have therefore been studying 
the resolved extra-planar stellar populatio ns of NGC 891 
|Mouhcine et al J 120071 ; iReikuba et alj|2008h . This galaxy is 
particularly interesting because (i) it is the nearest galaxy 
with similar morphology and mass to the Galaxy; (ii) previ- 
ous studies have presented no indications of disturbance in 
the stellar compone nt; and (iii) it is almost perfect ly edge-on 
(i = 89.8° ± 0.5° — K regel fc van der Kruiij .2005). essential 
to disentangle the halo and disk in the absence of kinematic 
information. 

One feature of NGC 891 that sets it apart from 
the Milky Way is that it possesses a substant ial H I 
halo (ISancisi fc Alien! 119791 ; lOosterloo et"al] 120071 ). which 
is due primarily to galactic fountain phenomena, although 
~ 10 % may be of accretion origin (|Fraternali fc Binnevl 
2008). It appears that this gaseous component is linked 
to active on-going star-formation, and indeed the ac- 
cretion rate matches closely the star -formation rate of 
2.9M Q yr _1 l|Fraternali fc Binnevl 120081 ). An ionised gas 
"halo" component, almost certainly related to the star- 



The distribution of stars in the Milky Way analogue NGC 891 3 



formation activity, has been detected as an Ha enve- 
lope surrounding the disk of the galaxy (|Dettmarl Il990l ; 
iRand et all Il990l : iRossa fc DettmaJ 12003 1. Star-formation 
in the Milky Way is known to be episodic, with nu- 
merous short-term peaks in the last 2 Gyr during which 
the star-formation rate reached 3-5 times the ste ady rate 
l|de la Fuente Marcos fc de la Fuente Marcos! 12004 ). Given 
the supernova-driven galactic fountain, it is likely that 
NGC 891 is currently undergoing such an episode. Since 
there is no reason to expect a causal connection between 
this halo gas and the halo or thick disk stellar populations, 
which in the Milky Way and in simulations are exclusively 
composed of ancient stars, we believe that NGC 891 can be 
examined as a plausible analogue of a galaxy formed in a 
similar manner to the Milky Way. 

As descr ibed in detail in a companion paper to this 
contribution (jReikuba et al.ll2008l . hereafter Paper II) , three 
fields in NGC 891 were observed with the Wide-Field Cam- 
era of the Advanced Camera for Surveys (ACS) on board the 
Hubble Space Telescope in cycle 11 (programme GO-9414). 
These three archival fields probe both the disk and thick disk 
of the galaxy over a wide radial range as well as the inner 
halo (see Fig. 1), and have slight overlaps which are very 
convenient as a means to ensure photometric consistency, 
and to check the reliability of the photometric uncertain- 
ties. With exposures of 7712 sec in both F606W (similar to 
broad-band V) and F814W (similar to broad-band 7) the 
photometry reaches to a limiting magnitude of / ~ 29, ap- 
proximately 3 magnitudes below the tip of the red giant 
branch (RGB). Paper II gives a thorough description of the 
data and the data reduction process, describing how the 
images were combined, and how crowded-field point-source 
photometry was measured from the image stacks; we there- 
fore refer the reader t o Paper II for th e se det ails. 

As we showed in iMouhcine et al.l |2007l . hereafter Pa- 
per I), the properties of the extra-planar stars in NGC 891 
questions the classical view of what a stellar halo should 
be. Our first analysis of the three contiguous ACS pointings 
showed that the metallicity distribution at ~ 10 kpc perpen- 
dicular to the disk of the galaxy is dominated by stars with 
[Fe/H] ~ —0.8. This is a full ~ 0.5 dex more metal-rich than 
what is measured for the Galaxy halo stars at a comparable 
height, and suggests that the metal-poor stars observed in 
the Galaxy may not be the domi nant population in galactic 
halos in general (IMo uhcine 20061). 

The purpose of the current contribution is to examine 
the spatial distribution of the stellar populations in more de- 
tail than in the previous two contributions in this series, and 
in particular to investigate the dumpiness of the halo popu- 
lation. The layout of this paper is as follows: in §2 we begin 
our analysis by reviewing the large-scale spatial structure as 
derived from 2MASS K-h&nd data; §3 presents the spatial 
distribution of the ACS star-count data, focussing on the 
minor axis profile and the vertical profile of the outskirts of 
the disk. In §4 we make a first attempt at uncovering small- 
scale sub-structures in the galaxy, and then proceed in §5 
to analyse the metallicity distribution over the survey area, 
and discuss the small-scale variations of the metallicity. The 
implications of these findings are discussed in §6, and we 
draw the major conclusions from our study in §7. 

Throughout this paper we assume a distance of 
9.73 Mpc (equivalent to a distance modulus of (m — M)o = 




Figure 2. The 2MASS if-band image of NGC 891 (upper panel), 
is fit within the area delimited by the mask (panel 'b') with a two- 
component disk plus bulge model (panel 'c'). The corresponding 
residuals are displayed in the bottom panel. 



29.94, from Paper I). We note that NGC 891 is located at 
relatively low Galactic latitude (£ = 140.38°, b = -17.42°), 
and therefore these fields suffer from significant (though not 
large) extinction fro m foreground dust: E(B — V) = 0.065 
l|SchIegel et all 19 98). corresponding approximately to Ay = 
0.22, A! = 0.13. 



2 GLOBAL MODEL FROM A'-BAND 
SURFACE PHOTOMETRY 

To obtain an insight into the global large-scale structure of 
the galaxy, we first perform a fit to the 2MASS near-infrared 
data on this galaxy. The upper panel of Fig. 2 shows the K - 
band 2MASS Large Galaxy Atlas data (Ijarrett et al.ll2003l ). 
retrieved from the NASA/IPAC archive, rotated 68° East 
of North to aid interpretation. In this Jf-band image the 
obscuring dust layer, which impedes easy analysis of the 
visib le DSS image i n Fig . 1, is very much attenuated. 

IXilouris et al.l (jl998l ) have previously made a detailed 
fit of this galaxy to A"-band data (and other wavelengths), 
using a composite model consisting of a disk, a bulge and a 
dust-layer. For the purposes of our own study, we felt that 
their model could be slightly improved on. In particular, 
they did not try to fit the boxiness of the NGC 891 bulge, 
which is very marked as one can immediately perceive in 
Fig. 2. For reference, their K-h&nA thin disk model has scale 
height hz = 0.32 ± 0.01 kpc and a scale length of Kr — 
3.93 ± 0.1 kpc, inclined at 9 = 89.6° ± 0.1° to the line of 
sight. 

We fitted the galaxy lig ht distribution usi ng the GAL- 
FIT software package v3.0 i|Peng et al.ll2002l ). The Large 
Galaxy Atlas images have 1" pixels (0.047 kpc), oversam- 
pling the original 2MASS frames by a factor of 2 x 2. For 
the purposes of determining t he noise characteris tics of the 
frame, we assumed (following iMaller et al.| [2005). that the 
image was derived from an average of 6 frames, each with 
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Figure 3. The upper panel shows a colour mosaic of all three ACS fields, which highlights the crowding and extinction problems close 
to the galactic plane, as well as the diffraction spikes, bright star halos and reflected light artefacts. The spatial distribution of point-like 
sources is displayed in lower panel, superimposed on the _R"-band 2MASS data. The sample has been culled by keeping only reliable RGB 
candidates, as detailed in the text. The selection region for the minor axis sample (Fig. 5) is marked in red (±10° from the minor axis), 
that for the outer disk (Fig. 7a) is marked green (16 < X < 18 kpc), while that of the "Solar-radius-like" sample (Fig. 7b) is marked 
blue (7 < X < 9 kpc). 



52 electrons read noise, and gain of 8.2 electrons per ADU. 
The disk was assumed to be a perfectly edge-on exponential 
(and hence Ki(R/hn) * (R/hn) in projection, where K\ is 
the modified Bessel- function) , with the following free param- 
eters: central surface brightness, scale length, scale height, 
position of centre and position angle. The bulge, however, 
was taken to be a Sersic model, with central surface bright- 
ness, effective radius r e , Sersic exponent n, axis ratio q (flat- 
tening), position of centre, and boxiness parameter c (see 
iPeng et al.ll2002l for a definition of this parameter). 

We limit the region to be fitted to \X\ < 15 kpc, as this 
corresponds approximately to the region where the disk fol- 



lows an exponential law (|van der Kruit fc Searldl 19811 ). and 
limit the fit in the vertical direction to \Z\ < 2.35 kpc (50") 
to avoid including regions obviously devoid of galaxy sig- 
nal in the 2MASS image (see Fig. 2). We masked out the 
stellar point sources in the frame (defined as 3-sigma devi- 
ations over the local background integrated over a 6.6 pixel 
aperture), blocking out circular regions of diameter three 
times the stellar FWHM. We also masked out the regions 
where the extinction is particularly strong and clumpy in 
the sub-area \X\ < 10 kpc, \Z\ < 0.55 kpc, which we identify 
as those pixels with H — K > 0.4. The dark pixels in panel 
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Figure 4. Panel 'a' shows the same data as Fig. 3, but with 
binned-up star-counts. The area under each of these 0.5 kpc X 
0.5 kpc pixels is calculated from the ACS image weight maps. 
We have additionally masked out the haloes of bright stars with 
appropriately-sized circular regions. Ellipses of minor axis length 
5, 7 and 9 kpc and axis ratio 0.73 are shown to demonstrate 
the presence of an extended elliptical component. The middle 
panel displays a de Vaucoulcurs model fit to the halo region at 
Z < —6 kpc, whose residuals are shown in the bottom panel. 



'b' of Fig. 2 show the areas that were masked out in this 
way. 

Curiously , when we test the K-band structural model of 
IXilouris et"afl (|l998l ) on the Jf-band 2MASS data, including 
their dusty disk model, but assuming that the galaxy is per- 
fectly edge-on, we do not find a particularly good fit. This 
slight variant of their model over-subtracts the thin disk, 
leaving behind a substantial excess of light at \Z\ > 0.25 kpc. 
In principle, this cannot be due to the minor difference in 
choice of disk inclination : the slight offs e t of 0. 4° between the 
edge-on model and the IXilouris et al.l J 19981 ) model would 
amount to a mere ~ 0.14 kpc apparent offset from Z = 
for stars in the galactic plane located at 21 kpc (the trunca- 
tion radius where the expone ntial profile drops precipitously 
Ivan der Kruit fc Searldll98ll ) either in front or behind the 
galaxy. However, the exponential model they use (and that 
we use too) is infinite in extent, leading to a substantial con- 
tribution at high \Z\ even for very small disk inclinations; 
this is not physical however, as the r eal disk is truncated at 
21 kpc (|van der Kruit fc SearleHl98lD . 

We therefore interpret the r e sidua l excess light at \Z\ > 
0.25 kpc from the lXilouris et all (|l998l ) model with 9 = 90° 



as a (slight) failing of that model due to the initial assump- 
tion that the disk is infinite. In the following we will also 
assume an infinite disk, but one that is perfectly edge-on; 
the outer regions are less problematic in this configuration, 
since their light is swamped by the much brighter inner re- 
gions of the disk. We feel that this assumption is justifie d 
by the V- and B-bands model fits of IXilouris et all l|l998h . 
which are consistent with the galaxy being perfectly edge- 
on. 



F or the reasons explained above, the IXilouris et al.l 
(1998) model probably underestimates the disk scaleheight. 
However, since their fitted value of hz = 0.32 ± 0.1 kpc is 
close to that of old stars in the Milky Way, it is informative 
to try to refit the galaxy with an additional thick-disk com- 
ponent to account for the excess light away from the galactic 
plane. Indeed, adding an additional disk to their model with 
a scale length constrained to be identical to their thin disk 
value (hit = 3.93 kpc) gives an excellent fit with a resulting 
scale height of hz = 0.76 ± 0.01 kpc. 

However, the data do not require two stellar disks; fit- 
ting with a single disk and bulge model gives an excellent 
reduced \ 2 (= 1-04, a value calculated assuming only Pois- 
son noise in the star-counts) . The fitted disk scale length and 
scale height are found to be ha = 4.19 ± 0.01 kpc and hz = 
0.57 ± 0.01 kpc, respectively. For comparison, t he Galactic 
thin disk has hn = 2.6 kpc and hz = 0.3 kpc (|juric et al.l 
2008). Note, however, that we have masked-out most of the 
disk close to the galactic plane (Fig. 2b), which must bias 
our scale height to large values. We found the best-fitting 
parameters to be: Sersic index n = 4.33 ±0.10 (close to a de 
Vaucouleurs profile), effective radius r e = 2.5± 0.1 kpc, flat- 
tening q — 0.73 ± 0.01, and boxiness c = 2.1 ± 0.1. Panel 'c' 
of Fig. 2 shows the resulting final model, and the residuals 
are displayed in the bottom panel. 

Although the residuals between the data and the model 
are quite small, there are large-scale patterns in the residual. 
These include the obvious thin structure close to the galac- 
tic plane, seen between —10 kpc < X < 10 kpc. Since this 
structure contains numerous point sources, which are clearly 
bright-star forming regions, it is likely a thin star-forming 
ring. The fact that it does not line up with the Z = axis 
and seems to undulate suggests an extra-planar influence in 
its formation. It is conceivable that this ring-like structure 
was formed in the interaction (or collision) with UGC 1809, 
recently claimed to be responsible for an extra-planar fila - 
ment of H I gas (|Oosterloo et al.ll2007l ; iMapelli et ai1l200ct ). 
Note that this extra-planar population lies on the opposite 
side of the Z = pla ne to what would be expected from the 
IXilouris et al.l l|l998h inclination angle of 9 — 89.6°, so it is 
not an artefact of our choice of assuming that the galaxy is 
perfectly edge-on. Another interesting aspect of the galaxy 
that is highlighted by the model fit in Fig. 2 is that the model 
over-subtracts the data beyond approximately 15 kpc; this 
corresponds to the break radius, where the observed pro- 
file b egins to drop much faster th an the inner exponential 
disk (Ivan der Kruit fc Searlelll98lf) . The fit also confirms the 
"boxy" appearance of the bulge, clearly seen in the upper 
panel of Fig. 2. 
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3 DISTRIBUTION OF POINT-SOURCES 

As can be seen in Fig. 1, there is significant overlap between 
the three ACS fields. Although we could have used the in- 
creased depth in the overlap regions to probe slightly deeper, 
we chose not to do so in order to maintain uniform sensitiv- 
ity over the survey region. Furthermore, to avoid problems 
with variations of source completeness, we chose to keep 
in the overlap areas only those data from the higher num- 
bered field (i.e. we kept the H3 data in the H2-H3 overlap, 
and the H2 data in the H1-H2 overlap). The bottom panel 
of Fig. 3 shows the distribution of point-sources detected 
with the D OLPHOT softw are package (a modified version 
of HSTphot, iDolphiiJliOQOh that have magnitudes I > 25.8, 
sharpness parameter less than 0.5 in both bands and crowd- 
ing parameter less than 0.3. This choice of sharpness and 
crowding parameters is identical to that adopted in Paper I. 
To be consistent with the colour-magnitude selection in Pa- 
per I (see Fig. 3 of that contribution), we impose a faint limit 
I < 27 A (for V -I < 1.3), then a linear increase to I < 27.3 
(at (V — I = 2.0) and finally a linear increase to I < 26.8 
(at V — I — 2.6). As we show in Paper II, this limit provides 
better than 50% completeness even in the relatively crowded 
regions at Z ~ — 4kpc. The spatial distribution of the 87973 
sources selected in this way is displayed in Fig. 3, superim- 
posed on the 2MASS A"-band image previously shown in 
Fig. 2. 

The irregular upper edge to the ACS point-sources in 
Fig. 3 (especially in fields HI and H2) is due to the imposed 
crowding limit; beyond this area the accuracy of the pho- 
tometry is compromised due to image blending. Also visible 
are the gaps between the ACS chips, as well as holes of low 
source density in the "haloes" around bright stars. 



3.1 The minor axis profiles 

A property of the stellar distribution that is immediately 
striking is the presence of a spheroidal distribution sur- 
rounding the bulge. This spheroidal population is seen more 
clearly if we bin the starcounts, as in Fig. 4, where we have 
overlaid ellipses of axis ratio 0.73 to guide the eye. The 
middle ellipse with minor axis length 7kpc delimits quite 
well the "edge" of this component. We stress that this is, of 
course, not a physical limit but simply an artefact of this 
particular greyscale representation. 

Is this ellipsoidal component the outer reaches of the 
bulge, or the inner parts of the halo? To answer this ques- 
tion we show in Fig. 5 the minor axis star-counts and the 
V and I-band surface brightness profiles. The star-counts 
profile (black histogram) was constructed by summing stars 
within a ±10° triangular wedge around the minor axis (be- 
tween the red lines in Fig. 3), accounting for lost detector 
area due to the various defects of the survey listed above. 
The blue and red lines in Fig. 5 show the surface bright- 
ness profiles from integrated light, measured in the F606W 
and F814W images respectively. These two profiles corre- 
spond to the median pixel value also in the same ±10° 
wedge, with a background value calculated from the interval 
—9 > Z > — llkpc subtracted off (this is to correct for the 
uniform "sky" background). These profiles are only plotted 
out to Z = —7.5 kpc, since at this distance the brightness 
is only ~ 1% of the background value; beyond this distance 
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Figure 5. The histogram shows the minor axis density distri- 
bution of the ACS star-counts, together with a fitted de Vau- 
couleurs model (green dashed line) which has effective radius of 
r e = 1.55 ± 0.03. The black line that is almost indistinguishable 
from the de Vaucouleurs profile for Z < —2 kpc is the minor 
axis profile of the galaxy model previously fit in Fig. 2. The blue 
and red lines show, respectively, the minor axis surface brightness 
profiles derived directly from the HST/ACS images. In all cases 
a ±10° wedge along the minor axis was selected for analysis (as 
indicated in the bottom panel of Fig. 3). 

the integrated light profiles are too uncertain. Finally, the 
black profile allows a comparison to the composite bulge 
plus disk model fitted in Fig. 2, again from a ±10° minor 
axis wedge. Since the normalisation between star-counts and 
surface brightness is unknown, we have shifted vertically the 
integrated light profiles, as well as the galaxy model profile, 
to agree with the star-count profile at Z = —5 kpc. 

It is patently clear from Fig. 5 that the ACS minor 
axis star-counts begin to become incomplete at Z > —4 kpc, 
and get progressively worse as one approaches the galactic 
planeQ Consequently, we ignore the histogram below \Z\ = 
4 kpc, and fit the profile of the ellipsoidal component with 
a de Vaucouleurs model in the range —4 > Z > —12 kpc; 
this yields an effective radius of r e = 1-55 ± 0.03 (the fit is 
shown with a green dashed line) The small difference be- 

1 This was the reason why we chose Z = —5 kpc as the position to 
normalise the surface brightness profile: further in, the star-counts 
are incomplete; further out, the uncertainties on the star-count 
distribution get larger. 

2 However, over this relatively small range in radius the profile 



The distribution of stars in the Milky Way analogue NGC 891 7 




-1 -2 -3 -4 -5 -6 -7 -8 

Z (kpc) 



Figure 6. The black line shows the distribution of (V-I) colour 
difference from integrated light, derived from a ±10° wedge on 
the minor axis. Given that the blue dip near Z ~ —5.5 kpc is 
likely an artefact due to scattered light, the profile is remarkably 
flat for Z < —2.5 kpc, indicating that the stellar populations and 
extinction are very uniform at large distances from the galactic 
plane. Correcting for the redde ning due to dust in the galaxy, as 
estimated from the H I map of lOosterloo et aL 1 1|2007|) , gives the 
red-line histogram. The error bars in both cases show the rms 
scatter in the median colour of the ACS pixels. 



tween this de Vaucouleurs fit to the "halo" star-counts and 
the disk plus bulge model (black line) fit to the global inner 
galaxy A'-band light is startling, and suggests an intimate 
connection between the halo and bulge. Note that the dis- 
crepancy between the two curves at Z > —2 kpc is due to 
the disk model beginning to contribute close to the galactic 
plane. 

The integrated light profiles are also very interesting. As 
one can see from Figs. 1 and 3, field HI just misses the cen- 
tre of the galaxy; this is the reason that the profiles begin at 
Z = —0.6 kpc. Evidently, in the region —4 > Z > —7.5 kpc 
both the V and I-band light follow each other and the star- 
counts profile extremely closely. This provides strong evi- 
dence that the stellar populations over this region are very 
similar, and that there is little variation in extinction. We 
probe this issue further in Fig. 6, where we show, with the 
black histogram, the distribution of colour difference be- 
tween the V and I-band profiles. Starting at Z — —7.5 kpc 
and proceeding towards the galactic plane, the colour dif- 
ferences are extremely small (consistent with zero) until 
Z = —5.5 kpc, where we encounter a large blue dip in the 
profile (centred at Z ~ —4.5 kpc), which clearly cannot be 
due to extinction. Inspection of the ACS images shows that 
this feature is likely an artefact due to scattered light (it 
is a diffuse elongated structure with a relatively sharp out- 
line and no obvious stellar counterpart). Beyond this dip, 



data can also be fit by an exponential function with scale height 
hz = 2.31 ± .04 kpc. Real galaxies can con tain significant halo 
sub-structure <MalirJll97sl:llbata et alj|2007l . e.g.), which renders 
the discussion of the statistical acceptability of different model 
profiles very uncertain over small radial ranges. 



the profile resumes at its previous level, but then reddens 
rapidly towards the galactic plane for Z > —2.5 kpc. 

Is this reddening of the colour profile towards the galac- 
tic plane due to stellar populations differences, or due to 
extinction? NGC 8 91 has significant amounts of dust out 
to \Z \ ~ 1.5-2 kpc (jHowk fc Savage! 1 19971 : iKamphuis et all 
l2007fl . so internal extinction appears to be a likely cul- 
prit. To assess this p ossibility, we use the H I maps of 
lOosterloo et al.l (|2007T h and assume a ratio of total hydro- 
gen column density to E(B — V) reddening of (5.94±0.37) x 
10 21 a toms cm" 2 mag -1 , id entical to that found in the Milky 
Way (|Rachford et al.ll2008l ). Correcting for half of this red- 
dening (i.e. assuming that half of the dust is in front of 
the stars, and half behind) gives the red-lined histogram in 
Fig. 6, which is consistent with a flat colour profile in the in- 
terval — 1 > Z > —3 kpc. We conclude, therefore, that there 
is no evidence from the colour profile for a change in stel- 
lar populations along the minor axis of NGC 891 between 
-1 > Z > -7.5 kpc. 

3.2 The thick disk 

The other striking feature of Figs. 3 and 4 is the presence 
of a thick envelope of stars surrounding the disk, reaching 
down from the Galactic plane to approximately Z ~ —4 kpc. 
This structure is particularly noticeable in field H3 between 
15 < X < 19 kpc. To explore this structure further, we 
show in panel 'a' of Fig. 7 the vertical star-counts pro- 
file at X = 17 ± 1 kpc with a black histogram. The blue 
and red curves show, respectively, the integrated light pro- 
files in the V and I-bands; as before we have truncated 
these at the point where they reach 1% of the background. 
The good correspondence between the integrated light pro- 
files and the star-counts profile is reassuring, though they 
highlight the incompleteness problems due to crowding at 
\Z\ < 0.5 kpc. Between —0.5 > Z > —5.0 kpc the star counts 
profile appears exponential; and fitting this region with an 
exponential model (green dashed line) gives a scale height 
hz = 1.44 ± 0.03 kpc. This value is very similar to the value 
of 0.9 ±0.18 kpc d educed for the thick disk of the Milky Way 
|juric et al.ll2008T ). and may be precisely such a component. 
This also will be discussed further in Paper II. The thin 
black line shows the profile of the galaxy model fitted to the 
2MASS data, and shown previously in Fig. 2, using the same 
normalisation as in Fig 5. Since the model fit to the 2MASS 
data does not incorporate a radial cut-off for the thin disk, it 
vastly over-predicts the counts in the inner ~ 1 kpc, and the 
absence of a thick disk component causes it to under-predict 
greatly the counts between —1 > Z > —5 kpc. However, be- 
yond Z ~ —5 kpc the model predicts a flat star-count distri- 
bution, rather close to what is observed. The flatness of the 
model is partly due to the fact that this vertical cut probes 
only a relatively limited extent in radius, but also partly 
because the model is a boxy spheroid. Indeed it is rather 
remarkable how well that bulge model, fitted primarily to 
the inner galaxy, works at this great distance. 

Panel 'b' of Fig. 7 shows an identical analysis for X = 
8±1 kpc, a location analogous to the Solar radius. Here again 
one can clearly see the presence of an exponential profile be- 
tween — 2 > Z > —5 kpc sandwiched between the thin disk 
and the halo. The green dashed line shows an exponential 
fit to the profile in this radial range, constrained to have a 
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Figure 7. Panel 'a': as Fig. 5, but for a vertical cut at X = 17 ± 1 kpc (the region contained between the green lines in Fig. 3). However, 
here a fitted exponential model with scale height hz = 1.44 ±0.03 kpc is shown instead. The black line shows the two-component galaxy 
model derived from 2-MASS data. As we discuss in the text, the model fails at Z > —5 kpc because it possesses a thin disk and no thick 
disk; it is astonishing however, that at Z < —5 kpc the 2MASS bulge model gives a close approximation to the outer halo star-counts. 
Panel 'b': identical to panel 'a', except that the vertical cut is made at X = 8±1 kpc. The green dashed line is an exponential, constrained 
to have a scale height hz = 1.44 kpc, and fit to the starcounts data in the range —2 > Z > —5 kpc. Clearly the "thick disk" component 
identified in the outer disk is present also in this region. 



scale height of hz = 1.44 kpc (the value measured from the 
data in panel 'a'); this shows nicely that the vertical profile 
of the thick disk is almost identical between X = 8 kpc and 
X = 17 kpc. 

Fig. 7 shows that we can extract a relatively clean "thick 
disk" sample by selecting stars between —2 > Z > —4 kpc: 
the upper limit is more than 3.5 thin disk scale heights above 
the galactic plane, while the lower limit lies before the halo- 
dominated region. The star-counts profile along two hori- 
zontal bands in this region is displayed in Fig. 8; the fitted 
projected exponential model was found to have a radial scale 
length of h R = 4.81±0.11 kpc in the -2 > Z > -3kpc band 
(upper panel) and ha = 4.83 ±0.15 in the —3 > Z > —4 kpc 
band (lower panel). These values are in excellent agreement 
with each other and are also similar to the radial scale length 
value of hit = 4.19 ± 0.01 kpc derived for the thin disk from 
the 2MASS Jf-band image in §2. These values should be re- 
fined in future work by performing two-dimensional fits to 
the disk region. However, this will likely require additional 
data near the galaxy plane (such data are already present 
in the archives, but not included in the present study since 
they are slightly shallower). 



3.3 2-D fit to outer halo 

Armed with the information of the spatial locations where 
the survey is reliably uncrowded and free of additional com- 
ponents, we can now attempt to fit a two-dimensional func- 
tion to the star-counts distribution in the halo. To this end, 
we use again the GALFIT software, fitting only the outer 
regions of the survey below Z — —6 kpc. Although we at- 
tempted to fit a general Sersic profile to these starcounts 
data, we found that the survey does not cover a sufficiently 
large radial range to fit both the Sersic index and scale 
length simultaneously. We therefore decided to fix the Sersic 
index to a value of n = 4 (i.e. a de Vaucouleurs profile). The 
resulting best fit (shown in the middle panel of Fig. 4) has 
scale length r e = 1.77±0.15 kpc, flattening of q = 0.50±0.03 
and boxiness of c = —0.46 ± 0.08. It should be noted that 
the fit is not good (reduced \ 2 = 1-41 with 449 degrees of 
freedom, which has a vanishingly small probability of occur- 
ring by chance). The residuals of the fit are displayed on 
the bottom panel of Fig. 4, and possess considerable non- 
symmetric large-scale features. Evidently the galaxy halo is 
not entirely symmetric. As we will see below in §4 and §6.4, 
there also appears to be small-scale spatial clumps in the 
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Figure 8. Horizontal cuts along the survey, showing the radial 
fall-off of the thick disk. The top and bottom panels correspond, 
respectively, to cuts between —2 > Z > — 3kpc and —3 > Z > 
—4 kpc. A projected exponential model is fit to the data in each 
panel (dashed line), resulting in a best-fit radial scale length of 
h R = 4.81 ±0.11 kpc and h R = 4.83 ±0.15 kpc for the cuts in the 
top and bottom panels, respectively. The gaps in the data (e.g. 
at 12.75 kpc in panel 'b') correspond to bins whose area is more 
than 20% missing or maskcd-out. 
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survey. It is these sub-structures (and to a lesser degree, the 
large-scale asymmetry), that are responsible for the poor \ 2 
of the fit. 



4 SPATIAL SUB-STRUCTURE 

Having studied the large-scale distribution of stars in 
NGC 891, we attempt next to search for localised clumps 
in the number density map of Fig. 4. The necessary first 
step is to subtract out the smoothly-varying components. 
We experimented using a two-dimensional Legendre polyno- 
mial fit (as performed on the metallicity profile in §5), but 
it transpires that the density gradient is too steep to give 
a satisfactory fit with low-order functions. We therefore re- 
turned to the GALFIT software to attempt a 2-dimensional 
fit of the star-counts distribution. 

To eliminate the effects of incompleteness due to crowd- 
ing, we masked the data in those regions with a stellar den- 
sity higher than 1000 stars kpc -2 (as suggested by the fits 
in Figs. 5 and 7). Since the star-counts distribution drops 
below this value (due to crowding) near the galactic plane, 



Figure 9. Panel 'a' displays a two-component (disk and spheroid) 
model fit to the star-counts distribution of Fig. 4; the corre- 
sponding residuals are shown in panel 'b'. As discussed in the 
text, the non-truncated disk model cannot fit well the region at 
X ~ 18 kpc, \Z\ < 3 kpc, which results in large scale residuals 
in this region. A ripple in the residuals centred near X ~ 5 kpc, 
Z ~ —7 kpc is also apparent. Panel 'c' displays the corresponding 
significance map, revealing several isolated pixels containing sig- 
nificant counts above the model; the same data is reproduced in 
panel 'd', where we have summed over the large Voronoi super- 
pixels used later for the metallicity analysis in Fig. 11. 

we had to rely on a model to define the boundary in sur- 
face density of the mask. However, since the fitted model 
depends on the boundary, we had to iterate a few times to 
achieve convergence. A two component galaxy with thick 
disk and halo was fit to the star-counts. This time, to allevi- 
ate the degeneracy of the problem, we fix the centre of both 
components to X, Z = (0, 0), we also fix the scale-length of 
the thick disk to ha — 4.83 kpc (as found in Fig. 8) and 
scale-height hz = 1.44kpc (as deduced from Fig. 7a), we 
further fix the halo to have a de Vaucouleurs profile and we 
set its boxiness/diskiness parameter to c = (thus the only 



10 R. Ibata, M. Mouhcine, M. Rejkuba 




Figure 10. The candidate satellites or satellite debris identified 
via their local deviation from the smooth model shown in Fig. 9. 
The middle and right-hand panels show the 0.5 kpc X 0.5 kpc pixel 
where the excess population is detected, while the left-hand pan- 
els display the colour-magnitude distribution of the stars in the 
region. To guide the eye we have overlaid three Galactic globu- 
lar cluster fiducial sequences, shifted to the distance modulus of 
NGC 891. These are, from left to right, M 15 ([Fe/H] = -2.2), 
NGC 1851 ([Fe/H] = -1.2) and 47 Tuc ([Fe/H] = -0.7). 

free parameters are the central surface brightness of the disk 
and halo, as well as the halo effective radius r e and flatten- 
ing q). We stress here that we are not attempting to fit a 
realistic galaxy model to the star-counts data, but simply 
trying to obtain a smooth representation of the data. (For 
completeness, we note that the fitted model had a halo ef- 
fective radius of r e = 1.97 ± 0.09 kpc, and a flattening of 
q — 0.55 ± 0.01, which are slightly different values to those 
fitted previously in §3.3). 

The residuals between the star-counts distribution 
(shown previously in Fig. 4) and this fit are displayed in the 
upper panel of Fig. 9. This residual map (panel 'b') shows 
large-scale structure that highlights the inadequacy of the 



simple model, fn particular, one can see a strong feature at 
X ~ 16 kpc, Z ~ —1.5 kpc which is due to the modelled 
disk not being sufficiently dense at that location. The fit 
cannot accommodate the absence of these stars slightly fur- 
ther out in radial distance (at X ~ 19 kpc, Z ~ —1.5 kpc) 
so compromises between the two. ft is possible that this is 
an artefact of fitting a model with a single exponential com- 
ponent. A further problematic region occurs at (X ~ 5 kpc, 
Z ~ —7 kpc), where one sees a large negative-positive ripple 
in the residuals (the positive part of the ripple being closer 
to the plane, at for instance, X ~ 4 kpc, Z ~ —5 kpc); curi- 
ously, the negative ripple is not reproduced on the symmetric 
side of the galaxy (X ~ —5 kpc, Z ~ —7 kpc). 

The significance of the residuals is displayed in the bot- 
tom panel, where we have assumed that the only source 
of uncertainty is the Poisson noise in the star-counts (i.e. 
that there is no noise in the fitting procedure). Although 
there are significant peaks around the previously-mentioned 
large-scale residuals, we also detect several significant iso- 
lated spikes. These isolated spikes are candidate satellites 
of NGC 891 or indeed accretion remnants. In Table 1 we 
list the positions of the 6 most significant (1-pixel) spikes 
above 2.5a that are well-separated from the areas with 
problematic residuals; Fig. 10 shows close-up views of the 
0.5 kpc x 0.5 kpc super-pixels in which they are found, to- 
gether with the corresponding colour-magnitude diagrams. 
Upon inspecting the ACS images visually, candidate satel- 
lite 'a' is found to possess a plausibly-concentrated dis- 
tribution (slightly off-centred of the super-pixel centre by 
SX r 0.05, SZ < 0.1), and a plausibly tight RGB se- 
quence. The nature of the other structures is less clear, al- 
though in all cases an RGB sequence is detected. 

The so-called "matched filter" method, where one uses 
a template stellar population (e.g., an RGB track) and a 
template spatial distribution (e.g., a 2-dimensional Gaus- 
sian) to optimally filter out contaminants, provides a more 
pow erful means to se arch for such spatial clumps (see, 
e.g., ISegall et al1l2007l . and references therein). The anal- 
ysis is somewhat complicated in the present situation by 
the strongly varying "background" population due to the 
smooth component of NGC 891, and also due to the complex 
spatial window function of the ACS survey. Furthermore, the 
sensitivity varies strongly as a function of the extent of the 
structure: for instance, one convincing satellite (not listed in 
Table 1) was detected visually, yet because of its compact- 
ness, and hence high stellar crowding, none of its RGB stars 
are present in the point-source catalogue. We have decided 
to defer a thorough matched-filter analysis to a subsequent 
contribution, since quantifying the detection and exclusion 
limits for the presence of sub-structure will require a lengthy 
discussion beyond the scope of the present contribution. 



5 SPATIAL DISTRIBUTION OF THE MEDIAN 
METALLICITY 

We now proceed to map out the metallicity properties of 
the surveyed stars. The metallicity of each star passing 
the quality tests laid out in §3 (with the additional con- 
straints I < 27.0, and DOLPHOT uncertainties SI < 0.1 
and 5V < 0.15), was derived by int erpolating between the 
models of IVandenBerg et all (|2006h of a-enhanced red gi- 
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Figure 11. Panel 'a': map of the median metallicity over the ACS survey area, colour-coded with the values shown on the wedge. A 
large-scale gradient in the halo region is noticeable, as well as a more metal-rich "thick disk". Panel 'b': map of the uncertainties on the 
median metallicity, derived from simulating the effect of the photometric measurement uncertainties on the metallicity interpolation. 
Panels 'c' and 'd' show the same information as 'a' and 'b', but with the maximum correction for internal extinction assuming that all 
the dust associated with the H I llOosterloo et al.ll2007l) is in front of the stars. 




Figure 12. The inhomogeneities in the median metallicity of the survey stars for the "standard" metallicities (left-hand panels) and 
with the maximum correction for internal extinction (right-hand panels) . The upper panels show a Legendre polynomial fit to the median 
metallicity map. The middle panels display the difference between the median metallicity map and this fit, while the lower panels report 
the significance level of the differences. These maps show that the halo has highly significant local variations in the median metallicity. 
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Table 1. Candidate sub-structures. 



Label in Fig. 10 X(kpc) Z(kpc) Significance (sigmas) 



a 8.25 -11.25 3.6 

b -2.75 -9.75 3.1 

c 0.75 -7.25 3.1 

d -4.75 -6.75 3.0 

e 11.25 -4.75 2.5 

f 16.25 -3.25 3.1 



ant branch stars of mass 0.8 M©, as discussed in Paper I. 
The models span -2.314 < [Fe/H] < -0.397; all stars out- 
side of the colour-magnitude range covered by these RGB 
tracks were flagged and not used in subsequent analysis (i.e., 
no extrapolation beyond the validity of the models was at- 
tempted) . 

After initially binning the survey area into 0.5 kpc x 
0.5 kpc pixels, we employed t he Vo ronoi tessellation algo- 
rithm of ICappellari fc Copinl (|2003f ) to implement a set of 
super-pixels in which the Poisson noise level of the star- 
counts exceeded S/N — 10. As we have discussed in Pa- 
per II, the RGB stars in all the galactic components iden- 
tified in NGC 891 cover a wide range in colour, implying a 
vast metallicity range. This intrinsic spread of stellar popu- 
lations makes the differences from population to population 
rather difficult to highlight, so we decided to make use of 
the median metallicity as our population statistic, as this is 
one of the most robust measures that can be made of a dis- 
tribution. This has the additional advantage that it is not 
affected in any way by the spatial density of the popula- 
tion, neither is it affected by issues such as inter-CCD gaps, 
holes due to bright stars etc, so in many ways it is one of 
the most sensitive and robust measures that one can make 
regarding the spatial properties of the galaxy. The map of 
the median value of the metallicity calculated in the Voronoi 
super-pixels is displayed in panel 'a' of Fig. 11. 

There are two main sources of uncertainty in calculat- 
ing the median metallicity map: a measurement uncertainty 
arising from the photon noise in the magnitude measure- 
ments of the stars, and a sampling uncertainty due to the 
limited number of stars per bin. Ideally, we would resample 
from the true metallicity distribution in each pixel, using 
the estimated photometric uncertainties to quantify the ef- 
fects of photon noise. While this would give the correct total 
uncertainty, it is not possible to implement, since the true 
[Fe/H] distribution in each pixel is unknown. 

We therefore proceed as follows to estimate the errors. 
To this end, we adopt the photometric measurement uncer- 
tainties returned by DOLPHOT; these estimates are reli- 
able, as we show in Paper II by comparing the photometry 



of stars detected in the overlapping region of fields H2 and 
H3.0 

• Step 1: The photometry of each of the Ni stars in super- 
pixel i is randomly drawn from a Gaussian distribution of 
dispersion equal to the measurement uncertainty (and with 
a mean equal to the actual measurement ), and its metallic- 
ity es timated by interpolation from the IVandenBerg et all 
(2006) tracks in exactly the same manner as for the sur- 
vey stars. This procedure is repeated 1000 times for each 
super-pixel. The resulting dispersion in the median values 
calculated from these Monte-Carlo simulations cri iroa i, using 
random realisations of the real stars actually present in each 
bin, gives a good estimate of the effect of the photometric 
measurement errors, but the sampling uncertainties have not 
yet been properly taken into account (the tilde is meant to 
highlight this fact). 

As mentioned above, we have no means of knowing the un- 
derlying [Fe/H] distribution in each super-pixel. However, 
we can calculate how much the Monte-Carlo procedure in 
"step 1" underestimates the total uncertainty, if the under- 
lying distribution were the observed distribution over the 
whole survey. 

• Step 2: To assess this, we draw an initial set of Ni stars 
at random from the full data-set, and then, using the corre- 
sponding photometric uncertainties, draw 1000 random re- 
alisations from this initial set, in an identical way to "step 
1". The dispersion of the median metallicity value cr^fake, 
is calculated from this sample. The whole procedure is re- 
peated 1000 times, in order to derive the average dispersion 
of the median (o^fake), that would be expected from apply- 
ing the procedure outlined in "step 1" to a super-pixel with 
Ni stars if there were no spatial stellar population gradients. 

• Step 3: Finally, we draw Ni stars at random from the 
full data-set (again with the photometry drawn randomly 
from a Gaussian distribution of dispersion equal to the corre- 
sponding photometric measurement uncertainties), and de- 
rive the median metallicity. The dispersion of the derived 
value over 1000 simulations gives the typical total uncer- 
tainty cr^fake in the median [Fe/H] in super-pixel i with Ni 
stars, which again is only strictly valid if the metallicity dis- 
tribution does not depend on spatial position. 

The difference between the result of "step 3" and "step 
2" gives the correction that needs to be applied to "step 1" 
to account properly for the sampling noise. Thus the final 
uncertainty on the median metallicity of bin i is calculated 
as Oi — a^roai + 0-^fakc — (?<,fake), where the last two terms 
are the corrections to the procedure outlined in "step 1". 
We note that these terms are corrections to an uncertainty, 
and therefore are not added in quadrature. The estimated 
uncertainties cr; are displayed in panel 'b' of Fig. 11. 

As discussed in §3, there is evidence for substantial 
amounts of dust close to the plane of the galaxy in NGC 891. 

3 We have also undertaken all of the following analysis adopting 
the alternative choice of using the mean magnitude-uncertainty 
relations derived in Paper II. The difference between the two 
methods of estimating the uncertainty is small and the con- 
clusions of this study remain identical. We chose to adopt the 
DOLPHOT error, because we judged that it was more reliable to 
have a star-by-star uncertainty estimate, especially for the outer 
halo region where crowding incompleteness is not an issue. 
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Figure 13. Panel a: the distribution of the statistical significance 
of the deviations from the model in panel 'c' of Fig. 12, showing 
a large number of highly significant super-pixels. The sample has 
a vanishingly small probability of being drawn by chance from a 
Gaussian distribution of unit dispersion. The filled-in histogram 
shows the pixels in the pure halo region below Z = — 5 kpc. Panel 
'b' shows the same information for the sample that was corrected 
for internal extinction (panel 'f of Fig. 12). This again cannot 
occur by chance. The dashed lines show a fitted Gaussian of dis- 
persion a = 1.4 (1.3) in the left-hand panel and a = 1.5 (1.4) 
in the right-hand panel, where the values in brackets refer to the 
halo-only sample. 



In an attempt to assess a nd correct the eff e cts of the dust, we 
again use the H I map of Oos terloo et al.l (|2007f ) to estimate 
th e E(B — V) reddening following the Galactic calibration 
of iRachford et al. (2008). In the following analysis, we will 
present the metallicity distributions with and without this 
additional extinction correction. Note that this is actually 
a maximum correction, assuming that all the extinction lies 
in front of the stars. The corresponding internal extinction- 
corrected median metallicity map and its uncertainty map 
are shown in panels 'c' and 'd' of Fig. 11. Although the cor- 
rection is substantial near the galactic plane (as can be seen 
in the corrected colour profile of Fig. 6) , the correction de- 
clines rapidly away from the plane, so that at the outer edge 
of the ACS survey region, the estimated internal extinction 
as derived in this way from the H I column density amounts 
to less than 0.001 mag. 

Fig. 11 (both panels 'a' and 'c') reveal a large-scale ra- 
dial gradient in the median metallicity of the halo, consis- 
tent with the results presented in Paper I; this can be seen 
clearly by inspecting the gradient along the X-direction at 
Z ~ —10 kpc. We perceive also that the median metallicity 
of the outer "thick disk" component is significantly more 
metal-rich than that of the halo, consistent with expecta- 
tions. We delay further discussion of the behaviour of the 
metallicity gradient to §6. 

The issue we wish to draw attention to here is the fact 
that the uncertainties on the median metallicities in the 
super-pixels are small, typically ~ 0.04 dex, compared to the 
large pixel to pixel variations. This can be seen explicitly in 
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Figure 14. Map of the AGB to RGB star ratio. The AGB star 
candidates are confined to regions close to the disk of the galaxy, 
indicating that the halo beyond \Z\ > 4 kpc shows no evidence 
for intermediate age stellar populations. 



Fig. 12; here the left-hand panels show the distributions not 
corrected for internal extinction, with the right-hand panels 
showing the internal extinction-corrected versions. Since we 
do not have an a-priori model of the large-scale distribution 
of metallicity, we adopt a pragmatic approach of fitting the 
observed distributions in Fig. 11 (panels 'a' and 'c') with 
a two-dimensional Legendre polynomial (up to order 3 in 
both X and Z, i.e. 10 parameters). The difference between 
the distribution of median metallicity and this "model" is 
displayed in the middle panels of Fig. 12, while the bottom 
panels show the significance level of these differences. An in- 
spection of the lower two panels reveals that the super-pixels 
have very significant localised variations in median metallic- 
ity that occur on the scale of a super-pixel or a few super- 
pixels. The actual distribution of the statistical significance 
of the variations in panels 'c' and T of Fig. 12 is displayed 
in Fig. 13, which demonstrates clearly that the variations 
are much larger than what one would expect from a halo 
with a smooth spatial distribution in median metallicity. 

We tested the statistical reliability of our analysis 
method by "scrambling" the photometry of the survey, i.e. 
keeping the positional information and super-pixel positions 
intact, but randomly reassigning the photometric measure- 
ments of the stars. This was undertaken only for those stars 
at Z < —5 kpc, which can be considered to represent a halo- 
only sample, that way we avoid mixing in more metal-rich 
populations from close to the galactic plane. The underlying 
smooth distribution was refit for each random realisation. 
Reassuringly, the resulting map of differences between data 
and model for this halo-only sample has a distribution of 
significance levels consistent with having been drawn from 
a Gaussian of unit dispersion with better than 10% proba- 
bility. 



6 DISCUSSION 

6.1 Median metallicity maps 

The median metallicity map (Fig. 11a), which shows both 
large-scale metallicity gradients and small-scale metallicity 
clumps, is one of the major results of this survey. However, 
before interpreting these maps as indicative of metallicity 
variations, it is important to consider the alternative: that 
they are revealing differences in mean age. Due to the well- 
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known age-metallicity degeneracy, a redder mean colour in a 
pixel can either correspond to a more metal-rich population 
or else an older population. One avenue to examine this pos- 
sibility is to examine the ratio of the spatial distribution of 
asymptotic giant branch (AGB) and RGB stars. Bright AGB 
stars tr ace predominantly intermed iate age (2 — 6Gyr) popu- 
lations l|Mouhcine fc Lanconll2002h , whereas RGB stars will 
predominantly represent more ancient populations. Their ra- 
tio therefore gives a relative map of the number of interme- 
diate to old stars in the galaxy. 

We selected candidate AGB stars to be those point- 
sources with 24.34 < Iq < 25.84 with all other qual- 
ity criteria identical to the RGB sample discussed in §4. 
As before, we amalgamated the 0.5 kpc x 0.5 kpc pixels 
in to super-pixels using the Voronoi tessellation algorithm 
of lCappellari fc Copinl (|2003l ). using a target signal-to-noise 
ratio of 10; the resulting map is shown in Fig. 14. There are 
substantially fewer AGB star candidates than RGB stars, so 
the algorithm is forced to generate very large super-pixels, 
especially in regions of low AGB density. The map shows im- 
mediately that the AGB star candidates are predominantly 
confined close to the disk of the galaxy. The completeness 
tests presented in Paper II show that at Z < —2 kpc the stel- 
lar distributions should not be significantly contaminated by 
blends. We therefore conclude that, since the intermediate 
age AGB stars appear not to be present in the halo regions, 
there is no evidence for variations in population age beyond 
Z < —4 kpc, and we therefore more confidently interpret 
colour variations as variations in metallicity. 

The colour scale in Fig. 11a was chosen to allow direct 
comparison t o the median meta llicity map for Galactic stars 
in the SDSS (jlvezic et al.ll2008l , their Fig. 18). This analysis 
of the SDSS covers much of the radial range 3 < R < 16 kpc, 
and vertical range 0.5 < \Z\ < 8 kpc in the Galaxy; and of 
course probes the Solar neighbourhood disk. In contrast, the 
present ACS star-counts survey is hopelessly incomplete in 
the thin disk, however it extends out much further in radius 
and vertical height. Essentially all of the metallicity gradi- 
ent in the Milky Way occurs at \Z\ < 4 kpc, which is only 
very poorly covered in the ACS fields HI and H2. Neverthe- 
less, we detect in fields H2 and H3 a clear vertical gradient 
between Z ~ —3 kpc and Z ~ —5 kpc, rather similar to 
what is seen in the Milky Way, albeit offset ~ 0.3 dex more 
metal-rich. 

We detect a significant drop in metallicity at large ra- 
dius, such that the outermost regions of field H3 are approx- 
imately ~ —0.2 dex more metal-poor than the inner halo at 
Z = —5 kpc on the minor axis (see also the discussion in Pa- 
per I). In contra st, the region of th e Galactic halo covered 
in the analysis of llvezic et al.l l|2008l ) is mostly uniform, with 
a slight positive metallicity gradient, due to the presence of 
the Monoceros ring structure at the upper radial extent of 
the surveyed region. 

The effect of the Monoceros Ring on the Galactic metal- 
licity distribution cautions against a simple interpretation of 
the differences in the halo components of the two galaxies. 
It is still possible that both NGC 891 and the Milky Way 
possess similar underlying ancient halos, that are radially 
well-mixed (over the scales surveyed) and metal-poor. By 
analogy to the Milky Way, the slightly more metal-rich in- 
ner halo (r 15 kpc) of NGC 891 may simply be due to 
a Monoceros Ring-like accretion event which is now better 



spatially mixed than the Monoceros Ring. However, if the 
stellar halo of NGC 891 is indeed dominated by intermediate 
metal-rich stars, the implications are profound: a moderately 
metal-rich stellar halo suggests an origin in the disruption of 
chemically evolved accreted fragments, which in turn has an 
implication on the mass function of small mass dark matter 
halos. Determining the dominant mechanisms that produce 
such a variation between galaxies of similar mass is a key 
prob l em for models of spiral galaxy formation l|Renda et al.l 
120051 ; iFont et al.ll2008h . 

The SDSS analysis presented by ICarollo et ail (I2007T ) 
suggests that the Galactic halo is composed of two distinct 
populations, an inner flattened halo with prograde rotation 
with a metallicity distribution that peaks at [Fe/H] ~ —1.5 
and an outer more spherical component, with retrograde 
rotation that peaks at [Fe/H] ~ —2.0 and that begins to 
dominate at radii between approximately 15 and 20 kpc. 
Inspection of the Legendre polynomial fits to the median 
metallicity in Fig. 12 (panels 'a' and 'd') reveals that this 
inner/outer halo dichotomy may be present in NGC 891 as 
well. In particular, the Fig. 12d, where we applied the cor- 
rection for internal extinction, shows that the halo metallic- 
ity is approximately constant everywhere except the outer 
edge of field H3 where the metallicity drops suddenly (at 
r ~ 15 kpc). Only a more spatially extensive survey will 
tell whether this reflects the global properties of a putative 
outer halo, or even whether the median metallicity continues 
to drop with distance to attain values consistent with those 
measured in the Galaxy. 



6.2 Metallicity clumps 

The analysis presented in §4 indicates that in addition to 
the large-scale variations discussed above, the median metal- 
licity varies significantly on the smallest scales that we can 
probe with the present dataset, ~ 1 kpc. This result is highly 
unexpected, especially for the halo component, since we ob- 
serve the galaxy in projection, and are integrating the pop- 
ulation over a large range along the line of sight. It is impor- 
tant therefore to consider all the possible explanations for 
this observation. 



6.2.1 Photometric zero-point? 

A first possibility is that the photometric zero-point is not 
constant over the ACS fields. This could be caused by, for 
instance, substantial flat-fielding errors. Simulations showed 
that, on average in the RGB selection region we employed, 
d[Fe/H]/d(V — /) = 1.08, so that very substantial variations 
in the zero-point (of order 0.1 mags) are required to in- 
duce the observed dumpiness in metallicity. The ACS is an 
extremely stable and well-calibrated instrument, and is rou- 
tinely used to derive exquisite photometry (5(V — I) < 0.01) 
of globular cluste rs which cover the field of view of the in - 
strument (see, e.g. lSaraiedini et al.ll2007l ; iPiotto et al.ll2007l ) . 
Clearly the ACS does not generally suffer from such flat- 
fielding errors. Furthermore, the observed metallicity clump- 
ing pattern in the present survey was not reproduced from 
field to field. This possibility can be eliminated. 
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Figure 15. The shift in reddening, compared to the adopted mean value of 0.065 from lSchlegel et al. 

I (Il998h . required to compensate for 
the apparent chemical inhomogeneities. Panel 'b' shows the same information for the sample that was corrected for internal extinction 



6.2.2 Photometric errors underestimated? 

A second possibility is that the photometric errors were 
poorly estimated. In this case the metallicity uncertain- 
ties would be correspondingly larger, diminishing the sig- 
nificance of the detections. We have run simulations to in- 
vestigate this possibility, and find that this can indeed cause 
the observed affect, though the uncertainties would have to 
have been underestimated by a factor of ~ 3. However, the 
photometry tests in the overlap regions presented in Paper II 
provide a very robust and reliable calibration of the uncer- 
tainties, and are consistent with the artificial star tests; we 
find it highly unlikely that the photometry errors were sig- 
nificantly underestimated, and a factor of 3 is completely 
out of the question. We judge therefore that this possibility 
can be also be safely eliminated. 



6.2.3 Foreground and background contamination? 

As we have discussed in Paper I, contamination from fore- 
ground or background sources is n egligible. Based on the 
Besancon Galaxy synthesis model (|Robin et al.l I2003T ). we 
expect ~ 119 stars in the range 24.5 < Iq < 27 in the entire 
survey area. The number of background sources is also very 
limited, as demonstrated by the Hubble Deep Fields North 
and South which contain, respectively, 10 and 22 point-like 
sou rces (FWHM < 0"2) w ith magnitudes 22 < F814W < 
26 (|Casertano et al.ll200d ). The ACS survey of NGC 891 
should therefore be essentially free of such contaminants. 
There should be, on average, less than one contaminant per 
Voronoi super-pixel in Fig. 11, and therefore foreground and 
background contaminants should have no significant effect 
on the median metallicity values. 



6.2.4 Foreground dust? 

It is also possible that we are detecting the effects of 
dumpiness in the foreground Galactic inter-stellar medium. 
Though probing other lines of sight, photometric studies of 
globular clusters with the Hubble Space Telescope, which 
show tight colour-magnitude sequences (especially at the 
main-sequence turnoff) again constrain the scatter around 
the mean extin ction. In the large sam ple of globular clus- 
ters studied by Saraicdi ni et al.l (|2007h . the typical scatter 
over an ACS field of view was o- E ( B _ V ) ~ 0.01 mag. In 



fields of higher extinction such as that of the globular cluster 
NGC 2808, where the mean reddening is E(B - V) = 0.23, 
accurate p hotometry constrains the differential re ddening to 
0.02 mag (|Bedin et al.l |2000| : iPiotto et al.ll2007r ). As men- 
tioned above, the mean reddening to wards NGC 891 is 
E(B - V) = 0.065 (|Schlegel et al.lll998h . so under the rea- 
sonable assumption that reddening variations scale approx- 
imately linearly with mean reddening, we expect less than 
one hundredth of a magnitude scatter in E(B — V) from fore- 
ground dust. 

In Fig. 15 we probe the shift in extinction, compared 
to the mean value of 0.065, that would be required to 
offset the deduced difference in metallicity. For this ex- 
ercise we simulated adding or subtracting extinction to 
each Voronoi super-pixel, re-computing the metallicity mea- 
surements on a star-by-star basis (extinction coefficients of 
A V /E{B - V) = 3.315 and Ai/E{B - V) = 1.940 were 
assumed). Although the rms dispersion of the map is only 
cr(E(B — V)) = 0.025 mag, it can be seen that some of the 
required offsets are huge, and are particularly problematic 
in metal-poor pixels, which occasionally would require im- 
possible negative extinction. For this reason, this pos s ibilit; 
also appears unlikely. Note that the ISchlegel et al.l (|l99i 
maps show relatively low variations in extinction towards 
NGC 891, as can be seen in Fig. 16. 



6.2.5 Dust in NGC 891 itself? 

NGC 891 is well-known to exhibit extra-planar dust fea- 
tures out to surprisingly large distances; indeed some of 
these structures can even be seen directly in the upper 
panel of Fig. 3. The vertical structures appear to be re- 
lated to galactic fountain or chimney phenomena, and driven 
by supernovae. Could dust in NGC 891 itself be caus- 
ing the observed variati ons in colour? Dust has been de- 
tected at \Z\ < 1.5 kpc (|Howk fc Savage! Il997r ). and molec- 
ular gas is also detected o ver a very similar vertical range 
|Garcia-Burillo et al . 1992). However, these distances are of 
course very close to the galactic plane in the context of the 
present study. 

In contrast, H I gas has recently been detected over es- 
sentially the entire ver tical range of the present ACS survey 
|Oosterloo et alJfeoOTh , forming a halo- like envelope around 
the galaxy. Earlier H I observat ions suggested that th e outer 
gaseous disk was highly flared (jSancisi fc Alledll979h or ex- 
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Figure 16. Map of the reddening in the vicinity of NGC 891 
(gnomonic projectio n centred on the galaxy), derived from the 
ISchleeel etail (1998) dust maps. The variations in the mean ex- 
tinction are clearly modest over this 1 sq. degree region, although 
the map obviously does not have the resolution required to probe 
the reddening variations over the ACS survey area (the marked 
polygon), much less to evaluate whether the small spatial scale 
variations suggested by Fig. 15 are realistic or not. 



tremely warped into the line of sight with the near side in 
the h emisphere probed by our survey l)Becauaert fc Combed 
1 19971 ). However, subseque nt modelling of the full velocity 
cube |Swaters et alj 1 19971 ) contradicted these conclusions, 
favouring a non-warped conf i gurati on. The very deep ob- 
servations of lOosterloo et all (|2007l ). which are among the 
deepest ever reached in any galaxy, show filaments and 
counter-rotat ing H I structures, sim ilar to those observed 
around M31 (|Braun fe Thilkerll2004h . A significant fraction 
of this gas is likel y to be of accret i on ori gin. 

The maps of lOosterloo et alJ l|2007h show that on the 
minor axis between Z — — 5kpc and Z — — 12kpc, the H I 
column density drops from 10 20 cm -2 to 10 19 cm -2 . If wc 
assume the same ratio of hydr ogen column density to ex- 
tinction as in the Milky Way (|Rachford et al.l l2008). this 
corresponds to a mean extinction of 0.017 to 0.0017 mag. 
However, given the contribution from accretions, the halo 
gas in NGC 891 is probably less abundant in dust than the 
Milky Way at the Solar neighbourhood, so the reddening in 
the halo of NGC 891 should be even lower than these esti- 
mates suggest. We therefore consider that dust in the target 
galaxy itself is unlikely to have contributed significantly to 
the observed small-scale variations in median colour beyond 
Z < -5kpc. 



6.3 Implications of the metallicity variations 

The remaining choice is then that the observed colour vari- 
ations imply real variations in metallicity. This is a fasci- 



nating possibility and would suggest that the ACS survey 
has revealed the existence of numerous accretion remnants 
that are spread over a large volume in the halo and are still 
far from being fully phase mixed. These remnants, which 
are identified by the change that they cause in the median 
metallicity of the average halo population, must nevertheless 
enhance the local density of the halo over the average value 
at the particular location of the structure. Why then (with 
the possible exception of the stream-like structure between 
X ~ lkpc, Z ~ — 9kpc and X ~ — 4kpc, Z ~ — 6kpc) is 
there no obvious correspondence between density variations 
(Fig. 9d) and metallicity variations (Fig. 11a)? This may 
be explained by the fact that it is much harder to identify 
density sub-structures in the presence of a strong density 
gradient (the surface density changes by more than an order 
of magnitude just in the halo region). The gradient itself 
will diminish the contrast of the structure, and errors in the 
model subtraction are likely to be much larger in a steep dis- 
tribution. In comparison, the large-scale metallicity gradient 
(Fig. 12a) is very gentle, which of course makes the survey 
considerably more sensitive for detecting localised variations 
in the median metallicity. 

The number of stars that are needed to cause the ap- 
parent metallicity variations turns out to be rather modest. 
Of course the actual fraction will depend on the metallic- 
ity distribution of the interloping population, but taking a 
somewhat extreme case of a single stellar population of uni- 
form metallicity [Fe/H] = -0.7 (like 47 Tuc), we find that 
in a typical Voronoi super-pixel a -0.1 dex metallicity vari- 
ation in Fig. 12b can be caused by a additional population 
numbering only ~ 15% of the putative "normal" smooth 
halo. If this is indeed the case, it is then not surprising that 
the population cannot be detected from its density contrast, 
since this is smaller than the change in density over a typical 
super-pixel. 

It is worth noting that due to the virtually negligi- 
ble contamination from foreground or background sources, 
the present observations are extremely sensitive to extended 
sub-structures. It is in principle possible with the ACS to 
detect a population of ~ 15 RGB stars scattered over a vol- 
ume of many hundred cubic kpc. Such sensitivity is orders 
of magnitude greater than what can currently be achieved 
in the Milky Way with the SDSS (although the SDSS per- 
forms extremely well at detecting concentrated structures, 
and has access to distance information). This would explain 
why sub-structures analogous to those apparent in Fig. 12 
are so hard to detect directly in the Galactic halo. 

The detection of such a highly (sub-)structured halo 
such as suggested by Figs. 12c and 12f has profound impli- 
cations for galaxy formation models. Hierarchical formation 
scenarios predict an abundance of su b-structures such as 
strea ms and merger remnants (see, e.g. Bull ock fc Johnston! 
2005). Indeed, the Via Lactea II simulation IjDiemand et al. 
20081 ) . the highest resolution simulation of the evolution 
of dark matter in a Milky Way like galaxy that has been 
published to date contains ~ 40000 dark matter satellites 
down to 10 6 Mq , some of which are tidally disrupted to 
form streams, especially in the inner regions of the modelled 
galaxy. It is possible that the chemical in-homogeneities re- 
vealed in Fig. 12 are caused by the disrupted remnants of 
this vast predicted population. 

Observationally, it is now well-established that the 
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Galactic halo contains numerous signatures of accretion 
( Bclo kurov et al. I l2006h . some of which, such as the giant 
stream of the Sa g ittariu s dwarf, are still on-going events 
l|lbata et al.lll994l . l2001bl lcl). Likewise, our giant neighbour 
Andromeda also d isplays a multitude of streams, shells, 
spurs and satellites ijlbata et al.ll2001al . l2007l : lFerguson et al.l 
2002). However, what is surprising about the results pre- 
sented here on NGC 891 is that the sub-structure appears 
ubiquitous: in the halo at Z < — 5kpc, 52% of the Voronoi 
super-pixels in Fig. 12 deviate by more than la from the 

smooth model. 

Although this appears astonishing, iBell et all (|2008l ) 
came to a similar conclusion in the Milky Way. These au- 
thors analysed the statistics of density deviations from a 
smooth model of the halo star-count distribution in the 
SDSS. They accomplished this by measuring the following 
statistic from a sample of halo stars: 

a/total = VnZ "' V 1 ; " 1 — , (i) 

n 2—ti 

where n is the number of discrete bins the sample is di- 
vided into, Di are the individual observed counts in those 
bins, Mi are the model values in the pixels, and M[ is a 
Poisson realisation of the model with mean Mi. Thus the 
nominator is the pixel scatter of the data around the model 
minus the expected scatter in the model, while the denom- 
inator was chosen to be the total of the Poisson scatter. 
This statistic is convenient because it is independent of the 
choice of bi nning scale, and the Poisson noise contribution 
is removed. IBell et all (|2008l ) found that that their statistic 
exceeds a/total = 0.4 in the Galactic halo region probed by 
the SDSS, and furthermore that projec ting the stream-only 
models of iBullock fc Johnston! l|2005l ) (i.e., those models 
with no underlying smooth halo component), and analysing 
them in the same way as the data, gave similarly large values 
of a/total. This indicates that the Galactic halo is substan- 
tially lumpy and also that it is approximately as lumpy as 
predicted by the most extreme hierarchical formation mod- 
els. 

The a/total statistic is related to the star-counts den- 
sity, and so not directly applicable to our median metal- 
licity maps. However, we can easily estimate the scatter 
that the a /total = 0.4 value corresponds to in a bin of 
100 stars (the target density for Voronoi super-pixels in the 
halo region in Figs. 11 and 12). Solving Eqn. 1 for Mi with 
Di = 100 and n = 1, gives model values of Mi = 59 or 
Mi = 142. We would therefore have interpreted the de- 
tection as a (Di — Mi)/\fD~i — 4.1a peak or equivalently 
a —4.2cr dip. This average scatter in density is substan- 
tially larger than the metallicity scatter that we report here, 
though this should be expected. First, our findings result 
from integrating along the line of sight through the entire 
halo, which must n ecessarily r e duce the contrast of any sub- 
structure, whereas IBell et al l l|2008h are able to break the 
halo into a (small) number of shells in distance. Second, the 
sub-structures will almost certainly have overlapping metal- 
licity distributions, which again must reduce the contrast. 
Thus we find that the present results on the variations of 
the median metallicity in NG C 891 are qualitatively similar 
to those of IBell et al.l (|2008l ) . and hence support hierarchi- 
cal halo formation scenarios, although further work is now 



needed to compare our observations to simulations, includ- 
ing the effect of chemical enrichment in halo sub-structures. 

It has not escaped our notice that the clumpy halo sug- 
gested by the metallicity variations we have just been dis- 
cussing appears inconsistent with the de Vaucouleurs profile 
that we have fit to the halo component of the galaxy in §2 
and §3. The de Vaucouleurs profile is known to be found in 
remnants of major mergers (see, e.g.. lNaab fc Truiilloll2006l . 
and references therein), and is associated with structurally 
well-mixed populations. While this may be the case for the 
smooth halo component, it is nevertheless possible that there 
are numerous additional small accretion structures that ar- 
rived after the main event that caused the formation of the 
bulk of the halo. We therefore judge that this apparent in- 
consistency is not too worrying, unless the majority of the 
halo is composed of sub-structure (but see the next section 
below) . 

6.4 Spatial dumpiness 

In the analysis in §4 and Fig. 9 we performed a simple 
model fit to the spatial distribution of RGB stars, which 
revealed certain systematic residuals as well as several candi- 
date dwarf satellite galaxies. The strong systematic residuals 
render useless most of the area of the map for the purpose of 
quantifying the incidence of low-level sub-structure, in the 
manner we have just undertaken in §6.3 with the distribution 
of median metallicity. However, the region at X > lOkpc 
and Z < — 4kpc appears devoid of nasty residuals or any 
very significant localised peaks in Fig. 9c. For this region 
we calculate a /tota l = 0.14, lower than the value obtained 
bv lBell et ail (|2008l ) for the Milky Way. The probability of 
finding such a large value of this statistic by chance is 0.8%, 
assuming the smooth model of Fig. 9a. As we have men- 
tioned above, we view a projection of NGC 891, whereas the 
IBell et alj (|200Sl ) resolved the Milky Way into various shells 
in distance. Thus one would expect a weaker sub-structure 
signal in NGC 891, consistent with what is seen. 

It should be noted that the ACS survey is very differ- 
ent to the SDSS, with ver y different photom etric problems. 
For instance, as far as the lBell et a l. (2008) analysis is con- 
cerned, one could worry about the photometric calibration 
of the SDSS over the sky, the effect of foreground and back- 
ground contamination, the consequences of the strategy of 
observing in long "stripes" , the effect of corrections for dust 
over large fractions of the sky, etc. The present ACS study 
has a different set of problems. The fact that we measure 
similar density fluctuations in both galaxies is therefore very 
reassuring for both analyses. 



6.5 Large-scale structures 

Despite the detection of a strong signal due to small-scale 
clumping, the spatial distribution of extra-planar stars in 
NGC 891 shows (Fig. 9) a surprisingly smooth large-scale 
distribution over a physical scale where the distribution 
of stars in M31 is highly clumpy with loops, spurs and 
streams of relati vely high surface brightness (nv ~ 27 - 
28mag/arcsec 2 ) |lbata et al.l 12001 j ; iFerguson et all 120021 ). 
It would appear therefore that the main stellar accretion 
events that NGC 891 has experienced were either not as 
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massive as those experienced by Andromeda, or that they 
happened longer ago and are now more smeared out. With 
the rather limited spatial coverage currently available for 
NGC 891, however, it is difficult to push this comparison 
much further. 



6.6 What is the relationship between the halo 
and bulge? 

In §2 we saw that fitting a two component model to the 
2MASS data at \Z\ ^ 0.5 kpc (the limit is approximate be- 
cause the mask in panel 'b' of Fig. 2 is irregular) gives 
a good fit with a bulge model that has a Sersic index 
of n — 4.33 ± 0.1, that is, a model that is essentially a 
de Vaucouleurs profile. Further out on the minor axis at 
Z < —4 kpc, a region which one would consider a-priori to 
be halo-dominated, both the integrated light profile and the 
star counts follow this same profile, and even the outer halo 
(§3.3) is consistent with the same r 1//4 -law. This begs the 
question whether there is a structural difference between 
the bulge and the halo. 

Unfortunately, the current data sets do not constrain 
the surface brightness pr ofile of the inner bul ge region sat- 
isfactorily. As argued bv lBalcells et al] (|2003l ). high spatial 
resolution data is required to measure accurately the Sersic 
index of bulges: seeing effects in ground-based surveys smear 
out point-sources in the nucleus of the galaxy, leading to er- 
roneous high Sersic index values. They find that bulges have 
Sersic indices n = 1.7 ± 0.7. Thus despite its relative immu- 
nity to extinction, the 2MASS T^-band data (with its ~ 2" 
resolution) is not particularly useful to probe the inner bulge 
profile. At the time of writing, the Hubble Space Telescope 
archive has a single high-resolution NICMOS/NIC3 image 
in F160W (i.e. approximately H-band) that covers the very 
centre of the galaxy. However, even this NICMOS image suf- 
fers significantly from the strong interstellar extinction that 
blankets the centre of NGC 891, rendering any study of the 
central light profile highly uncert a in. 

In light of the Balcells et al.l (|2003f ) study, the paucity 
of reliable data in the inner galaxy forces us to interpret 
our findings cautiously. The boxy bulge in NGC 891 seen 
clearly in Fig. 2, has most probably in reality a Sersic pro- 
file with a low (exponential-like) index, consistent with the 
modern picture of these structur es as having been f ormed 
from disk instabilities (see, e.g. I Athanassoulal [2008) . The 
fact that we derive de Vaucouleurs-like index values in the 
fit to the _K"-band 2MASS data is likely due to the fact that 
we have included data out to \Z\ < 2.35 where the halo 
component begins to dominate. The reasonable match of 
this same model to the outer halo profile (at Z < —5 kpc 
in Fig. 7a) indicates that this de Vaucouleurs-like profile is 
a reasonable model for the entire halo region probed in the 
ACS survey. The ellipsoidal structure detected in panel 'a' 
of Fig. 4, and fitted with a "disky" elliptical structure in 
§3.3, therefore need not have a direct relationship with the 
inner boxy bulge. Nevertheless, it will be interesting to ob- 
tain further high-resolution infra-red imaging of the centre 
of this galaxy to actually measure the Sersic index of the 
bulge. 



6.7 Thick disk 

Thick disks are now known to be a ubi quitous component 
of galaxies (|Yoachim fc Dalca nton 2006;|). These structures 
may hold important clues to the early formation of disk 
galaxies, as they may be in themselves fossil evidence of the 
formation mechanism. Currently it is not clear whether these 
structures form by secul ar processes that redistribute star s 
in the galaxy (see, e.g., iKroupal [20021 ; iRoskar et all 120081 ). 
whether they are the result of the heating of th e early thin 
disk via a significant merger (|Quinn et al.lll993T ). or else due 
to repeated ac cretions of low- mas s satellites over the lifetime 
of the galaxy (lAbadi et alj|2003h . 

NGC 891 is an excellent target to study the thick disk 
component: it is edge-on and close enough to resolve indi- 
vidual stars, and far enough to afford us a panoramic view. 
The thi ck disk of NGC 891 was first d etected in integrated 
lig ht by van der Kruit fc Searlei (|l98ll ) , and later confirmed 
by iMorrison et al.l } 19971 ). who found that this component 



possesses a scale height between 1.5-2.5 kpc, and comprises 
between 3% and 21% of t he surface brigh t ness o f the thin 
disk at the galactic plane. IMorrison et all (|l997f ) conclude 
that the thick and thin disks have approximately twice the 
scale height of the Milky Way disks, which seems to be 
consistent with the larger extent of the H I and molecu- 
lar gas in this galaxy perpendicular to the disk. The thick 
disk structure has also been detected from the mid-infrared 
emission of asymptotic giant branch stars, which p ossess a 
scale-height of 1.3 ± 0.3 kpc (|Burgdorf et al.ll2007l ). How- 
ever, the present contribution provides the first detection 
of the thick disk in this galaxy from star-counts, which 
gives much higher signal to noise for fitting purposes. It 
is interesting that the thick disk component is structurally 
very close to a perfect exponential disk, with a scale height 
hz ~ 1.44 ± 0.03 kpc at X = 17kpc which remains essen- 
tially constant to X = 8 kpc, and with a scale length of 
h R = 4.8 ±0.1 kpc, larger, but rather close to the value of 
4.19 ± 0.01 kpc deduced for the thin disk. For comparison, 
the Galactic thick disk has a scale length of h R = 3.6 kpc, 
substantially longer than the h R = 2.6 kpc scale length of 
the thin disk, implying a substantial scatter between galax- 
ies in the relative properties of the two disk components. 
The well-behaved exponential implies a well-settled compo- 
nent in dynamical equilibrium. We defer a discussion of the 
implications of this finding to Paper II, where we will addi- 
tionally investigate the stellar populations in the component 
as well as their profile. 

I n a sample of edge-on galaxies studied in integrated 
light, lYoachim fc Dalcanton! (2006) find the following scal- 
ing relationships between radial scale length and circular 
velocity V c : h R ,±hin = 3.4(V C /100 km s -1 ) 1 ' 2 kpc, h R<t hick = 
3.9(V C /100 kms -1 ) 10 kpc . Given the measured circular ve- 
locity of V c = 224 km s" 1 (|Paturel et al.ll2003l ). these rela- 
tions predict h R , t hin = 8.95 kpc and h R , t hick = 8.74, com- 
pletely inconsistent with our measurements in NGC 891. 
It is curious that these relations also do not work for the 
Milky Way; it will be interesting to investigate why these 



4 We have refrained from re-analysing the relative masses of the 
thick and thin disk, since new ACS data cover the disk at X ~ 
12 kpc, and we feel that this question can be addressed much 
better with that additional field. 
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two galaxies, which after all are among the best cases one 
can stud y, seem to deviate so ma r kedly from the global trend 
found bv lYoachim fc Dalcantonl (|2006h . 



7 CONCLUSIONS 

We have undertaken a structural and colour (i.e. chemical) 
analysis of a large survey of resolved stars in the edge- 
on galaxy NGC 891. Extending out to X ~ 24kpc and 
Z ~ — 12kpc, the survey provides a unique view of the stel- 
lar populations in the outer bulge, halo and thick disk of a 
galaxy whose optical properties are very similar to those of 
the Milky Way. 

Our main findings can be summarised as follows: 

(i) We detect the thick disk of NGC 891 in star-counts, 
and give the first measurement of the global structural prop- 
erties of the component based on these measurements, which 
are inherently much more reliable than earlier measurements 
based on integrated light. We find that in the spatial window 
of this survey, the thick disk follows closely an exponential 
disk profile with scale length of hjt — 4.8±0.1 kpc. The scale 
height is found to be hz = 1.44 ± 0.03 kpc. The smooth ex- 
ponential nature of the component suggests that it is now 
dynamically well mixed. We will discuss the stellar popula- 
tion properties of the component further in Paper II. 

(ii) The large-scale metallicity distribution is approxi- 
mately uniform in the halo component out to ~ 15 kpc, 
with a median value of [Fe/H] ~ —1.1. Furthermore, the 
colour of this inner region of the halo is identical to that 
of the bulge beyond ~ 0.5 kpc (after accounting for inter- 
nal reddening due to dust in NGC 891 itself). At a ra- 
dial distance of approximately 15 kpc the halo metallicity 
drops to [Fe/H] ~ — 1.3, and in analogy to the Milky Way 
jCarollo et al.1 120071 ). this is suggestive of the presence of 
a distinct outer halo component. Further data beyond the 
present survey limits are required to confirm this possibility. 

(iii) Along the minor axis, the structural profile of the 
halo merges smoothly into that of the bulge, with the com- 
bination of the two following closely a de Vaucouleurs profile 
from ~ 0.5 kpc out to the edge of the survey. However, the 
shape of the distribution clearly changes with distance, be- 
ing rounder (q — 0.73) close to the bulge, and substantially 
flatter (q = 0.50) at large distance. The similarity of the 
bulge and halo in both colour and structural properties is 
striking, though we reserve judgement on the relationship 
between the two due to the difficulty of studying the struc- 
ture of the bulge region inside 0.5 kpc. 

(iv) The most important finding of this work is the dis- 
covery of small-scale sub-structure in the median metallicity 
maps (Figs. 11 and 12). Though the metallicity variations 
are small, they are nevertheless highly significant, and ap- 
pear to indicate that the halo of NGC 891 is composed of 
a multitude of small accretions that have not yet been fully 
blended into a "normal" smooth halo. This result based on 
metalli city variat i ons, is in broad agreement with the find- 
ings of iBell et al.l l|2008l ). based on densit y variations in the 
Milky Way. We also repeat the analysis of IBell et all (120081 ) 
on the stellar density distribution in the halo, and find that 
their a /total statistic, the fractional rms deviation of the 
data from a smooth model, is cr/total = 0.14. This is some- 
what lower than the value of 0.4 found in the Milky Way, 



but may be consistent with it given that we view NGC 891 
in projection. 

(v) We detect directly 6 low-mass satellite candidates 
from their local density enhancements; these are possible 
dwarf spheroidal galaxies, or tidally disrupted fragments 
thereof. However, further work is needed to implement a 
more sophisticated detection algorithm (such as a matched 
filter), to quantify the detection limits and physical prop- 
erties of the satellites, and to compare these to hierarchical 
formation scenarios of low-mass structures in a A-CDM cos- 
mology. 

Thus we have found that the halo of NGC 891 shows 
evidence for substantial amounts of sub-structure. These 
clumps show similar statistics in the scatter of stellar den- 
sity as the Milky Way, a finding which supports the most 
extreme models of hierarchical halo formation. We have also 
detected the presence of these sub-structures from the small- 
scale statistical fluctuations in the metallicity distribution, 
the first-time it has been possible in any galaxy, and which 
gives complementary information on the nature of the sub- 
structures. Substantial follow-up work is now required to 
undertake cosmological galaxy formation simulations, incor- 
porating star-formation and chemical evolution in the dark 
matter mini-halos and to compare the resulting stellar dis- 
tributions to our observations. 

The current ACS survey is unfortunately rather limited 
in spatial extent and many of the questions left open by this 
study will require a more panoramic view to be addressed. 
It will be interesting also to extend this analysis to other 
extra-galactic systems, and in particular to assess the level 
of sub-structure as a function of galaxy morphology, type, 
mass, and accretion history. 
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